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We propose a theoretical scheme to generate a controllable and switchable coupling between two double- 
quantum-dot (DQD) spin qubits by using a transmission line resonator (TLR) as a bus system. We study dy- 
namical behaviors of quantum correlations described by entanglement correlation (EC) and discord correlation 
(DC) between two DQD spin qubits when the two spin qubits and the TLR are initially prepared in X-type 
quantum states and a coherent state, respectively. We demonstrate that in the EC death regions there exist DC 
stationary states in which the stable DC amplification or degradation can be generated during the dynamical 
evolution. It is shown that these DC stationary states can be controlled by initial-state parameters, the coupling, 
and detuning between qubits and the TLR. We reveal the full synchronization and anti-synchronization phenom- 
ena in the EC and DC time evolution, and show that the EC and DC synchronization and anti-synchronization 
depends on the initial-state parameters of the two DQD spin qubits. These results shed new light on dynamics 
of quantum correlations. 

PACS numbers: 03.67.Lx, 42.50.Ct, 03.67.Bg , 73.21.La 



I. INTRODUCTION 

It is well known that quantum entanglement 01 0] and 
quantum discord [3] are two different types of quantum cor- 
relations. In recent years, it has been widely recognized that 
both of the two quantum correlations are essential quantum 
resources which can be used to realize quantum information 
processing, and the fact that quantum discord is a more gen- 
eral concept to measure quantum correlations than quantum 
entanglement since there is a nonzero quantum discord in 
some separable mixed states 1141- 4 1 ill . However, the question 
of the relation between entanglement correlation (EC) and dis- 
cord correlation (DC) is still an open problem in the field of 
quantum correlations. 

Interactions of quantum systems are at the core of quan- 
tum information processing. In particular, quantum comput- 
ing requires that inter-qubit interactions are controllable, and 
can be selectively switched on and off Ill2lll3ll . Many phys- 
ical systems have been explored for the realization of practi- 
cal quantum information processors, such as cavity quantum 
electrodynamics (QED) system [14], optical system IU5ll . and 
solid-state system B16I422I1 . Solid-state devices are the promis- 
ing candidates for the implementation of quantum computa- 
tion due to the possibility of fabricating large integrated net- 
works. Among different kinds of solid systems, double quan- 
tum dot (DQD) system 123143 ill and a transmission line res- 
onator (TLR) 132l438ll are particularly attractive because of 
the relative long spin coherence time and high controllabil- 
ity of DQD system and quantum bus function of the TLR. 
Several proposals have been proposed to realize controllable 
couplings and local operations of quantum dot (QD) qubits via 
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a TLR 1I39I - I42I1 . In Ref. 63] , the electron spins in nanowire 
QDs couple to the electric component of the resonator elec- 
tromagnetic field and enable quantum information processing 
in an all-electrical fashion. In Ref. ll40ll . all-electrical cou- 
pling between QD spin qubits and a TLR are also used to 
produce effective interactions between spin qubits. In Ref. 
B4lll . two QD spin qubits were embedded in a superconduct- 
ing microstrip cavity, virtual photons in a common cavity 
mode could mediate coherent interactions between two distant 
qubits. In this paper, we want to propose a new scheme to im- 
plement the controllable coupling between double-quantum- 
dot (DQD) spin qubits in terms of the magnetic coupling be- 
tween DQDs and the TLR by using the TLR as the quantum 
bus system. In our scheme, each DQD spin qubit couples to 
the TLR through the magnetic filed generated by the current of 
the TLR. We will study the relation between EC and DC by in- 
vestigating the dynamic behaviors of EC and DC between the 
two DQD spin qubits in the combining system consisting of 
two DQDs and a TLR. We will show the appearance of the DC 
stationary states in the EC death regimes and demonstrate the 
full synchronization and anti-synchronization of DC and EC 
in the time evolution. In particular, it shall be indicated that 
the DC between the two DQD spin qubits can reach and keep 
its maximum fixed even if the EC disappears completely in 
the time evolution for certain initially prepared X-type states. 

This paper is organized as follows. In Sec. |IlJ we present 
our coupling scheme of two DQD spin qubits by using the 
TLR as the quantum bus system, and show that the inter-spin- 
qubit interaction in our coupling scheme is controllable and 
switchable. In Sec. [Till we investigate dynamical behaviors 
of EC and DC between two DQD qubits in the controllable 
coupling scheme, and discuss the relation between EC and 
DC for X-type initial states. The stationary amplification and 
degradation of the DC and the time synchronization and anti- 
synchronization of the EC and DC are revealed in the time 
evolution of the combined hybrid system. Finally, we con- 
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elude this work in Sec.lIVI 



II. CONTROLLABLE COUPLING SCHEME OF TWO 
DQD SPIN QUBITS 

In this section, we propose a scheme to generate an effec- 
tive controllable interaction between two DQD spin qubits by 
using a TLR as the data bus. The combined system under our 
consideration is indicated in Fig. 1. It consists of two DQDs 
charged with two excess electrons and a TLR. The length of 
the TLR is L. The two DQDs, 1 and 2, are placed at the po- 
sitions L/4 and 3L/4 of the TLR, respectively. These are the 
antinode of the quantized current of the TLR. Two electron 
spins in each DQD are localized in adjacent QDs, coupled via 
tunnelling. The distance between the lower dot and the TLR 
and the distance between the two dots are r. 

In terms of the annihilation and creation operators a and a ' , 
we can write the Hamiltonian for the TLR as 
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FIG. 1: The proposed setup with two DQDs, biased with external 
potential A^, / = 1,2, magnetically coupled to a TLR of length L. 
The two DQDs 1 and 2 are placed at the positions L/4 and 3L/4 
of the TLR, respectively. These are the antinode of the quantized 
current of the TLR. In each DQD, both of the distance between the 
lower dot and the TLR and the distance between the two dots are r. 
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where oj r is the frequency of the TLR. The Hamiltonian of a 
DQD is most conveniently written in the two-electron singlet- 



triplet basis {ISnXirO),!^), 
tion axis in the z-direction as 
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and the auxiliary singlet state with two electrons in one quan- 
tum dot, IS02), is coupled via tunneling t to the separated 
singlet, |5 n), where the subscript m, n denotes the dot occu- 
pancy. Es and Ej are the energy of the \S\\) and 17^) states, 
respectively. Ao is the energy difference between the \S n) and 
|5 02) states set by the electric field. B e is the external mag- 
netic field, fiB is the Bohr magneton and gs the electron spin 
g-factor. 

If the external magnetic field in the upper dot B" is different 
from that in the lower dot B[„ and if we consider the effect 
of the hyperfine interaction with the nuclear spins which can 
be studied by adding a static frozen effective nuclear field B U N 
(Z?' w ) at the upper (lower) dot to the total magnetic field 11431 - 
14511 . the interaction between the spin and the total magnetic 
field in each DQD is given by Hi - -gB(*B(B„ ■ S „ + B/ ■ S i)/H 
which can be rewritten as 

Hl = _^£ ( B i( + B /) . ( 5 i( + 5 ;)/2 , 
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where the total magnetic fields which an electron in each DQD 
experiences are B u — B U N + B" and B/ = B' N + B' e , the sum op- 
erators and the difference operators of the two electron spins 
can be expressed in terms of the two-electron singlet-triplet 
basis as 
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Eqs. and (O indicate that the homogeneous part of the 
magnetic field B u + B/ simply adds vectorially to the external 
field B e , changing slightly the Zeeman splitting and preferred 
spin orientation of the triplet states. The inhomogeneous part 
AB = B u — B[, on the other hand, couples the triplet states 
{|r?i>, 177,), IT")} to the singlet state |S n >. 

The degree of mixing between two states will depend 
strongly on the energy difference between them. In the case 
of gBfiBB e » t and gs^B V(AB 2 ), the spin-aligned states \T^ X ) 
and \T U ) are split off due to Zeeman energy gB/UBB e , then the 
perturbation of these states will be small, while the spin-anti- 
aligned state \T X1 ) remains mixed with the state |Sn). Then 
we can write the Hamiltonian of this combined system as 
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H = haj r a' ] ' a + hj^^&l + (A^ + E J S )\S 02 ) j(S Q2 \ 
7=1 

+ y B fi B (AB z j )di + t j (\Su)j(So2\ + \So2)j(Su\), (6) 
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where a>j = (E J S - E J T )/h, AB Z . is the z-component of ABj, and 
we have introduced Puali spin operators 



cri = \Sn)j{Su\-\T° n )j{T^\ 
cri = \Sn}j(T n \ + \T u }j(Sn\ 



(7) 



In Hamiltonian ^ we have used only 17^)^ and \Sn)j as 
a spin-qubit degree of freedom. Qubit-TLR Interaction is in- 
duced naturally by the magnetic field gradient AB Z which in- 
cludes the contributions from the nuclear magnetic field and 
the TLR itself 



A2K = A* + 



4nr ' 
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where ;Uo is the vacuum permeability, and AB Z N is the gradient 
of the longitudinal component of the nuclear magnetic fields 
between the two dots of the jth DQD. The current Ij at the 
position of the jth DQD due to the resonator of length L is 
quantized as 



(9) 



where I is the inductance per unit length of the TLR. 

Substituting Eqs. (O and (0 into Eq. (|6]l we arrive at the 
following Hamiltonian 
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In the interaction picture with respect to the term (Ai + 

E J S )\S 02) j(S 02], we obtain the interaction Hamiltonian after 
discarding rapidly oscillating terms 
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When ojj » gs^B^B z N ,/(2fi), and ojj + oj r » ojj - oj r , 

S %mr yHt' tne rotating-wave approximation can be applied 
to get an effective Hamiltonian of the combined system 
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H 2 = haj r a l 'a + tiY + (-ly^hgia&i + &ia f ), (12) 

7=1 

where the effective coupling constant g = sj^jj-- Eq. 

< TT~2T > is the Hamiltonian of the usual Jaynes-Cummings model 
of two atoms with &{ = |5n>/r°| and ai = |r°>/Sn|. 



If both DQDs are strongly detuned from the TLR, i.e., \6j\ - 
jj - u> r \ » \g\, we can adiabatically eliminate the TLR mode 



using the following transformation 



U = exp 
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To second order in the small parameters g/(5/, the effective 
Hamiltonian becomes 



H e f f = na» r fl a + — > 
7=1 
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where the last term describes swap of the DQD states through 
virtual interaction with the TLR. 

To easily switch on and off the inter-qubit coupling is one 
of the most important open problems in quantum computing 
hardware. Here we propose a way to overcome the severe 
problem. From Eqs. ([8]), ( fT2b and ( TBI ) we can see that the 
inter-spin-qubit coupling between the two DQDs is mediated 
by the TLR, and the interaction between each DQD and the 
TLR originated from the magnetic field gradient fioIj/4nr be- 
tween the two dots of each DQD, which is produced by the 
TLR. So if additionally we add an asymmetric magnetic field 
B. in the z-direction, then we have 



AB Z 



(15) 



where dBi is the difference of the asymmetric magnetic field 
B, between the two QDs of the jth DQD. Then by tuning the 
magnitude of dBi, we can get AB'j = 0. This implies that 
the interaction between the jth DQD and the TLR is turned 
off. Therefore, we can conclude that the effective coupling 
between the two DQDs in the present spin-qubit coupling 
scheme is controllable and switchable. 



m. DYNAMICS OF QUANTUM CORRELATIONS 
BETWEEN TWO DQD SPIN QUBITS 

In this section, we investigate dynamics of quantum cor- 
relations between two DQD spin qubits. We will study dy- 
namic evolution of quantum entanglement and quantum dis- 
cord when the two DQD spin qubits are initially in three- 
parameter two-qubit X-type quantum states [46] which play 
an important role in a number of physical systems. 

Firstly, we solve the Hamiltonian of the DQD-TLR system 
under our consideration. In the interaction picture with respect 
to the first term hoj r a^a of the Hamiltonian 77 e ff d!4l . we get 



< H = hY J Clj&l - hxicrld 2 + d-la-l). 



(16) 
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where the inter-qubit and qubit-TLR couplings are given by 
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with N - a* a being the number operator of the TLR. We can 
write «# in the basis {|S>i|S> 2 , |S>i|r> 2 , \T),\S) 2 , |7*>i|r>2} as 



( Q.y+Q.2 

Qi - Q 2 -x 
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-Q, -D. 
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Here and in the after we use \S)j and \T)j to replace ISn); 
and \T^)j in Sec. [TTJ respectively. The four eigenstates of the 

Hamiltonian < H can be obtained as 



l^l> = |S>l|S> 2 , 

PF 2 > = cos(0/2)|5>,|r> 2 + sin(0/2)|r> 1 |5> 2 , 

PF 3 > = -sin(e/2)|5) 1 |r> 2 + cos(0/2)|r>,|5) 2 , 

1^4) = |r>iir> 2 , 



(19) 



with the corresponding eigenvalues 

Ei = -E 4 = S(Qi + Q 2 ), 



E 2 = -E^h^-^+x 2 - 
The mixing angle in Eq. ( fT9] > is defined by 

sin = -nx/E 2 , cos = - Q 2 )/£ 2 . 



(20) 
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From Eqs. ( fT9l and (f20b we get the time evolution oper- 
ator ofHasU = exp(-/£ , „f/^)|T„)( , P, 1 | which can be 
expressed as 
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with the following parameters 
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We assume the two DQDs are initially prepared in a class 
of state with maximally mixed marginals (pA(B) = M(B)/2) de- 
scribed by the three-parameter X-type density matrix p(0) = 

(Iab + XlLi c '®"a ® &*b) 1^ wmcn can t> e expressed as 



P(0) = - 



' 1 + C3 Ci - c 2 

1 - C3 ci + c 2 

ci + c 2 1 — C3 

\c x - c 2 1 + c 3 J 



(24) 



where c,- (0 < |c,| < 1) are real numbers satisfying the unit 
trace and positivity conditions of the density operator p. The 
density operator p includes the Werner states and the Bell 
states as two special cases. Here and in the after we use the 
label A(B) to denote the first (second) DQD. 

If we initially prepare the TLR in the coherent state 



l<M0)> = \a) r 



(25) 



Then, the reduced density matrix of the two DQDs at time 
t can be obtained as 



pit) = Tr r (t/ T p(0)®|^(0)X<A r (0)|t/). 



(26) 



For simplicity but without loss generality, here we consider 
the case of two identical DQDs, i.e., we take u>a - <^>b — <*> 
and 6a = 8b = 8. Then from Eqs. ( TP71 ) and (f2Tb we have 
cos — 0. In this case, the explicit form of density operator at 
time t is 



Pit) = ~ 



(I + c 3 

1 — C3 C\ + c 2 

ci + c 2 1 - C3 
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where we have introduced the following parameter 

(ci - c 2 ) e 2i(u+ « 2/,5) ' exp [-|a| 2 (l - e 4i « 1,/s )]. 
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(27) 
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In what follows, we study dynamic properties of quantum 
correlations between the two DQD spin qubits in terms of 
the expression of the density operator of the two DQD spin 
qubits given in Eq. ( 1271 ). We investigate in detail EC and DC 
described by concurrence [47] and quantum discord, respec- 
tively. When the density matrix of the two qubit system has 
an X-type structure, the concurrence has a simple analytic ex- 
pression B46I1 



C(t) AB =2 max{0,Ai(0,A 2 (t)h 



(29) 



where K x {t) = \p u (t)\ - ^p 22 (t)p^(t) and A 2 (f) = |p 23 (f)| - 
^pn(t)pu(t) with p,j (i,j = 1,2,3,4) being the matrix ele- 
ments of the density operator pit) in Eq. (27). For the density 
operator given by Eq. (27) we have 



Ai(f) = -Gel - |1 - c 3 |), 
A 2 (0 = j(|ci+c 2 |-|l+c 3 |). 



(30) 



Now we turn to investigate dynamic evolution of quantum 
discord correlation between two DQD qubits under our con- 
sideration. Quantum discord [3] is defined as the difference 
between the total correlation and the classical correlation with 
the following expression 



D(p)=I(p A :p B )-C(p). 



(31) 



Here the total correlation in a bipartite quantum state p is mea- 
sured by quantum mutual information given by 



I (Pa ■ Pb) = S (p A ) + S (p B ) -S (p), 



(32) 



where S (p) — -Tr(plogp) is the von Neumann entropy, 
Pa = Tib(P) and ps = Tia(P) are the reduced density oper- 
ators for subsystems A and B, respectively. And the classical 
correlation between the two subsystems A and B can be de- 
fined as 



C(p) = max 



S(pa)-Y,P* S (p ( a) 



S(p A ) - min 
{Pt) 



(33) 
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Here {P/J is a set of projects performed locally on the sub- 
system B, and p* = ^Tr B \^t A ® p (l A ® P^j is the state 
of the subsystem A conditioned on the measurement outcome 
labeled by k ,where pt — Tr^ ]^I A ® Pk) P {J a ® Pkj\ denotes 
the probability relating to the outcome k, and I A denotes the 
identity operator for the subsystem A. 

In order to obtain the quantum discord for the two qubit- 
system [48], we first evaluate the mutual information of state 
p(t) given in Eq. ( fSTb . The four eigenvalues of p(t) are 



U,2 



-(i-c 3 ±h +c 2 |), 



^3,4 = ^(1 +C3 ±N). 

Then the mutual information reads 



I (Pa : Pb) = 2 + ^ log 2 A;. 



(34) 



(35) 
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Note that here we have used S (p A (t)) - S (pb(0) = T since the 
two reduced density matrixes p A {t) and ps(f) are maximally 
mixed, that is p A (t) = I A /2 and psif) = 

For calculation of the amount for the classical correlation 
C(p) defined in Eq. (33), we propose the complete set of or- 
thogonal projectors {Pf, = \6k)(0k\,k =11, -L} for a local mea- 
surement performed on the subsystem B, where the two pro- 
jectors are defined in terms of the following two orthogonal 
states 



|0|[> = cos^|0> + e^sin^|l), 
\6 ± ) = e^sin(/j|0)-cos^|l), 



(36) 



with < ip < n/2 and < (p < 2n. After the two project mea- 
surements with p|| = p x = 1 /2, the reduced density matrices 
of subsystem A read as 



„H _ 1 / 2(1 + C3 cos 2ip) esin2y> 

^ A ~ 4\e*sin2y> 2(1 - C3 cos2i/p) 



1 / 2(1 - C3 cos 2ip) -esin2^ 



^ A 4 1 -e* sin2i,c> 2(1 + C3 cos2i/p) 

where we have introduced the parameter 
e = (ci + c 2 )e-y + c e ! '< 



(37) 



(38) 



According to Eq. (37), it is straightforward to obtain the 
eigenvalues of the reduced density matrix p* as follows 



4 a = & a = 2 (1± ^> 



where we have defined T as 



r = -1/C3 cos 2 (2(/j) + '—j- sin z (2i^>), 



(39) 



(40) 



then we have 



S{p\) = s( PA -)=f(T) 



l+r, /i+r\ i-r, /i-r\ 



which leads to the classical correlation 

C{p{t)) = l-min[/(D]. 

v4 



(42) 



Since the function /(T) is a monotonically decreasing func- 
tion, in order to get the minimal value of /(F) we should 
choose proper parameters ip and (p to ensure the parameter F 
defined in Eq. (40) is maximal. From Eqs. (38) and (40) it is 
easy to get the following inequality 



r < 



yjc^ cos 2 (2tp) + c 2 A sin 2 (2(/?) 



C4, 



for IC3] > C4, 
for |c3| < C4. 



where we have introduced the parameter 

|ci + c 2 | + |c | 

C4 = — - — . 

If we define y(i) as 

y(t) = max[|c 3 |,c 4 ], 
then the classical correlation can be expressed as 



C(p(f)) = 



1 + (-i) m r 



io g2 [i + (-iy r ]. 
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Therefore, the quantum discord can be written as 
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D(p(t)) = 2 + J] ^ lo 82 ^' - C(p(0), 



(43) 



(44) 



(45) 



(46) 



(47) 



where the amount of the classical correlation C(p(f)) is given 
by Eq. (46). In principle, we have obtained the dynamics of 
the quantum discord according the above expression given in 
Eq. (47), provided that we know the initial condition of the 
system. In what follows we will study dynamic properties of 
the quantum entanglement and discord for some initial states 
in detail. 

From the density operator of the two DQD spin qubits at 
time t d27| i and Eq. (45) we can obtain the parameter 



r(0 



|c 3 |. 

C4, 



for 2|c 3 | > |ci - c 2 ] + |ci + c 2 \. 



for 2|c 3 | < |ci - c 2 |e 



-2|of 



+ |Cl +C 2 \. 



(48) 



which indicates that the classical correlation expressed by Eq. 
(46) largely depends on the initial state parameters of the two 
DQDs and the TLR {c1.C2.c3, a]. Then making use of Eqs. 
(46-48) we can obtain expressions of the quantum discord in 
different regimes of the initial state parameters. When 2|c3| > 
l c i ~ c i\ + ki + c i\, we nave 

4 2 
£>(p(f)) = 2 + ^log 2 /l,--2 2 5>g2r m , (49) 



!=1 

and when 2|c3| < \c\ - c 2 \e 
4 



-2|n|- 



+ \c\ + c 2 |, we have 
2 



D(p(t)) = 2 + ^i,log 2 ^-^^log 2r ;„, (50) 
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FIG. 2: (Color online) Time evolution of concurrence (dot dashed 
curves) and quantum discord (solid curves) for different values of C3 
when other parameters are taken by c\ = 1, c% = -C3, a = 2, g = 1, 
and 6=10. 



FIG. 3: (Color online) Time evolution of concurrence (dot dashed 
curve) and discord (solid curve) for different values of a. Other pa- 
rameters are c\ = 1, ci = -C3 = -0.3, g = 1, and 6 = 10. 



where two parameters y m and y' m are defined by 

7m = 1 + (-i)»|c 3 |, y m = 1 + (-iy C4 . (5i) 

In order to see clearly dynamic characteristics of the quan- 
tum entanglement and quantum discord, in the following we 
numerically investigate the time evolution of the concurrence 
given by Eqs. ( l29j > and d30l >. the discord D(p{t)) given by Eqs. 
(49) and (50). 

Fig. 2 indicates the influence of initial states of the two 
DQD spin qubits on dynamic evolution of the EC and DC 
when the DQD spin qubits and the TLR are initial X-type 
quantum states and a coherent state, respectively. We have 
plotted time evolution of concurrence (dot dashed curves) and 
quantum discord (solid curves) for some X-type initial quan- 
tum states of the two spin qubits when other parameters take 
fixed values. Fig. 2 reflects some interesting dynamic proper- 
ties of EC and DC. 

(i) Fig. 2 demonstrates that in the EC death regions there 
exist DC stationary states in which DC does not change in a 
finite time area. From Fig. 2(a-d) we can see that in the DC 
stationary states the stable DC amplification or degradation 
can be generated during the dynamical evolution. We can also 
see that these DC stationary states can be controlled by initial- 
state parameters, the coupling, and detuning between qubits 
and the TLR. Fig. 2(a-b) clearly indicate that the DC may be 
amplified in the most time of the death regime of the EC, and 
remains its maximal value. However, Fig. 2(c-d) reflect the 
fact that the DC may be degraded in the most time of the death 
regime of the EC, and remains its minimal value. These imply 
that one can generate stable DC amplification (degradation) in 
the EC death regions for certain initial X-type states. 

(ii) Fig. 2 indicates the appearance of the full synchroniza- 
tion and anti-synchronization phenomena in the EC and DC 
time evolution, and shows that the EC and DC synchroniza- 



tion (anti-synchronization) depends on the initial-state param- 
eters of the two DQD spin qubits. Fig. 2(c) and Fig. 2(d) 
correspond to the EC and DC synchronization evolution while 
Fig. 2(b) corresponds to the EC and DC anti-synchronization 
evolution. From Fig. 2 we can also see that both EC and DC 
time evolutions have the same period which is independent of 
C2 and C3. Indeed, when c\ — 1 and 1 > C3 = —ci > the 
evolution period of EC and DC can be analytically obtained 
from Eqs. (28-30) with the simple expression T = 6n/(2g). 
This implies that the larger period can available by increasing 
(decreasing) the detuning 6 (the effective coupling constant g). 

(iii) Fig. 2 indicates that both EC and DC may exhibit sud- 
den death phenomenon under certain conditions in their dy- 
namic evolution. Fig. 2(a) shows that the discord sudden 
death (DSD) periodically occurs while EC always remains 
zero. From Fig. 2(b) we can see that the entanglement sudden 
death (ESD) periodically happens while DC almost remains 
its maximal value in the EC death regions. 

In Fig. 3 we show the effect of initial states of the TLR 
on EC and DC dynamics. We have plotted time evolution of 
concurrence (dot dashed curves) and quantum discord (solid 
curves) for some initial quantum states of the TLR when 
other parameters take fixed values. From Fig. 3 we can see 
that the time evolutions of EC and DC are completely anti- 
synchronous with the same period, and the ESD occurs for 
enough large values of the initial parameter a in the EC dy- 
namics. From Fig. 3(c) and Fig. 3(d) we can also see that the 
increase of the initial-state parameter a of the TLR not only 
leads to the ESD but also the appearance of the DC station- 
ary states in the EC death regimes. The lifetime of these DC 
stationary states can be increased by increasing value of the 
initial-state parameter a. 

In Fig. 4 we investigate the influence of the detuning pa- 
rameter on dynamic evolution of the EC and DC. We have 
plotted time evolution of concurrence (dot dashed curves) and 



7 



o 
o 

C/5 



CD 
O 



CD 



3 
O 



O 

o 




FIG. 4: (Color online) Time evolution of concurrence (dot dashed 
curves) and quantum discord (solid curves) for different values of 6 
when other parameters are taken by c\ = \,Ci = -C3 = 0.3, a = 2, 
and g = 1. 



quantum discord (solid curves) for different values of the de- 
tuning of the TLR for a fixed initial state of the combined 
system and a fixed coupling. From Fig. 4 we can see that the 
lifetime of the DC stationary states and the EC death regimes 
increases with the increase of the detuning parameter 6. The 
larger the detuning, the longer the lifetime of the DC station- 
ary states and the EC death time become. By using Eqs. (28- 
30), similar numerical analysis indicate that the lifetime of the 
DC stationary states and the EC death regimes can be manipu- 
lated by changing the coupling constant g. Therefore, we can 
conclude that the lifetime of the DC stationary states and the 
EC death regimes can be controlled by changing the detuning 
parameter and the coupling constant between the DQD and 
the TLR. 



IV. CONCLUDING REMARKS 

In conclusion, we have proposed a magnetic coupling 
scheme between two DQD spin qubits by using the TLR as 
a bus system. We have shown that the inter-spin-qubit cou- 



pling is controllable and switchable. The coupling control- 
lability can be realized by changing external magnetic field. 
It is worthwhile to mention that the present scheme can be 
used to build a hybrid qubit system where DQD spin qubits 
are integrated together with superconductiong qubits in the 
same TLR. We have studied in some detail dynamical behav- 
iors of inter-qubit quantum correlations described by EC and 
DC when the two DQD spin qubits and the TLR are initially 
prepared in some X-type quantum states and a coherent state, 
respectively. We have studied the relation between EC and 
DC. We have demonstrated that in the EC death regions there 
exist DC stationary states. In these DC stationary states, the 
stable DC amplification or degradation can be generated dur- 
ing the dynamical evolution. We have shown that the lifetime 
of these DC stationary states depends on the initial-state pa- 
rameters, the coupling, and detuning between qubits and the 
TLR. We have also found the full synchronization and anti- 
synchronization phenomena in the time evolution of the EC 
and DC, and indicated that the time synchronization (anti- 
synchronization) of the EC and DC dynamics depend on the 
initial-state parameters of the two DQD spin qubits. We have 
indicated that both EC and DC may exhibit sudden death 
phenomenon under certain conditions in their dynamic evo- 
lution, and the DC stationary states always appear in the EC 
death regions. The DC stationary states, synchronization and 
anti-synchronization found in the present paper highlight new 
characteristics in the EC and DC dynamics of the combined 
DQD-TLR system, and have potential applications in quan- 
tum information processing. It could be convinced that the 
magnetic coupling scheme suggested in the present paper can 
be used to realize quantum information processing since it 
takes advantages of controllable and switchable coupling and 
qubit scalability. 



Acknowledgments 

This work was supported by the National Fundamental Re- 
search Program Grant No. 2007CB925204, the National Nat- 
ural Science Foundation under Grant Nos. 10775048 and 
1 1075050, the Program for Changjiang Scholars and Innova- 
tive Research Team in University under Grant No. IRT0964, 
and the Research Fund of Hunan Provincial Education De- 
partment Grant No. 08W012. 



[1] L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. 

Phys. 80, 517 (2008). 
[2] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, 

Rev. Mod. Phys. 81, 865 (2009). 
[3] H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 

(2001). 

[4] B. Schumacher and M. D. Westmoreland, Phys. Rev. A 74, 
042305 (2006). 

[5] B. Groisman, S. Popescu, and A. Winter, Phys. Rev. A 72 
032317 (2005) 



[6] L. Henderson and V. Vedral J. Phys. A: Math. Theor. 34, 6899 

(2001); V. Vedral , Phys. Rev. Lett. 90, 050401 (2003). 
[7] J. Oppenheim, M. Horodecki, P. Horodecki, and R. Horodecki, 

Phys. Rev. Lett. 89, 180402 (2002). 
[8] D. Yang, M. Horodecki, and Z.D. Wang, Phys. Rev. Lett. 

101,140501 (2008). 
[9] D. L. Zhou, Phys. Rev. Lett. 101, 180505 (2008). 
[10] D. Kaszlikowski, A. Sen(De), U. Sen, V. Vedral, and A. Winter, 

Phys. Rev. Lett. 101, 070502 (2008). 
[11] M. Piani, P. Horodecki and R. Horodecki , Phys. Rev. Lett. 100, 



8 



090502 (2008); M. Piani, M. Christandl, C. E. Mora and P. 
Horodecki , ibid. 102, 250503 (2009). 
[12] Y. Makhlin, G. Schoen, A. Shnirman, Nature (London) 398, 
305 (1999). 

[13] Y. X. Liu, L. F. Wei, J. S. Tsai, F. Nori, Phys. Rev. Lett. 96, 
067003 (2006). 

[14] C. J. Hood, M. S. Chapman, T. W. Lynn, and HJ. Kimble, Phys. 

Rev. Lett. 80, 4157 (1998) 
[15] P. Kok, W.J Munro, K. Nemoto, T.C. Ralph, J. P. Dowling, and 

G. J. Milburn, Rev. Mod. Phys. 79, 135 (2007). 
[16] Y. Makhlin, G. Schon and A. Shnirman, Rev. Mod. Phys. 73, 

357 (2001). 

[17] R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha and L. 

M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007). 
[18] M. Blencowe, Nature (London) 468, 44 (1010). 
[19] I. Chiorescu, N. Groll, S. Bertaina, T. Mori, and S. Myiashita, 

Phys. Rev. B 82, 024413 (2010). 
[20] D. I. Schuster, A. P. Sears, E. Ginossar, L. DiCarlo, L. Frunzio, 

J. J. L. Morton, H. Wu, G. A. D. Briggs, B. B. Buckley, D.D. 

Awschalom, and R. J. Schoelkopf, Phys. Rev. Lett. 105, 140501 

(2010). 

[21] Y. Kubo, F. R. Ong, P. Bertet, D. Vion, V. Jacques, D. Zheng, 
A. Dreau, J.-F. Roch, A. Auffeves, F. Jelezko, J. Wrachtrup, M. 
F. Barthe, P. Bergonzo, and D. Esteve, Phys. Rev. Lett. 105, 
140502 (2010). 

[22] H. Wu, R. E. George, J. H. Wesenberg, K. Molmer, D. I. Schus- 
ter, R. J. Schoelkopf, K. M. Itoh, A. Ardavan, J. J. L. Morton, 
and G. A. Briggs, Phys. Rev. Lett. 105, 140503 (2010). 

[23] V. N. Golovach, A. Khaetskii, and D. Loss, Phys. Rev. Lett. 93, 
016601 (2004). 

[24] A. V. Khaetskii and Y. V. Nazarov, Phys. Rev. B 61, 12639 
(2000). 

[25] H-A. Engel and D. Loss, Phys. Rev. B 65, 195321 (2002). 

[26] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, K. C. 
Nowack, T. Meunier, L. P. Kouwenhoven and L. M. K. Vander- 
sypen, Nature (London) 442,766 (2006). 

[27] K. C. Nowack, F. H. L. Koppens, Y. V. Nazarov, and L. M. K. 
Vandersypen, Science 318, 1430 (2007). 

[28] V. N. Golovach, M. Borhani and D. Loss, Phys. Rev. B 74, 
165319 (2006). 

[29] E. A. Laird, C. Barthel, E. I. Rashba , C. M. Marcus, M. P. 

Hanson, and A. C. Gossard, Phys. Rev. Lett. 99, 246601 (2007). 
[30] J. T. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, 

M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 

Science 309,2180 (2005). 
[31] T. Fujisawa, T. H. Oosterkamp, W. G. van der Wiel, B. W. 



Broer, R. Aguado, S. Tarucha, L. P. Kouwenhoven, Science 
282, 932 (1998); T. H. Oosterkamp, T. Fujisawa, W. G. van 
der Wiel, K. Ishibashi, R. V. Hijman, S. Tarucha, L.P Kouwen- 
hoven, Nature (London) 395, 873 (1998). 

[32] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin and, R. J. 
Schoelkopf, Phys. Rev. A 69, 062320 (2004). 

[33] A. Blais, J. Gambetta, A. Wallraff, D. I. Schuster, S. M. Girvin, 
M. H. Devoret, and R. J. Schoelkopf, Phys. Rev. A 75, 032329 

(2007) . 

[34] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.- S. Huang, 

J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Na- 

ture(London) 431, 162 (2004). 
[35] D. I. Schuster, A. A. Houck, J. A. Schreier, A. Wallraff, J. M. 

Gambetta, A. Blais, L. Frunzio, J. Majer, B. Johnson, and M. 

H. Devoret, Nature (London) 445, 515 (2007). 
[36] Q. Q. Wu, J. Q. Liao, and L. M. Kuang, Chin. Phys. Lett. 25, 

1179 (2008). 

[37] J. Majer, J. M. Chow, J. M. Gambetta, Jens Koch, B. R. John- 
son, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, and 
A. Wallraff, Nature (London) 449, 443 (2007). 

[38] M. A. Sillanpaa, J. I. Park, and R. W. Simmonds, Nature (Lon- 
don) 449, 438 (2007). 

[39] M. Trif, V. N. Golovach, and D. Loss, Phys. Rev. B 77, 045434 

(2008) . 

[40] Z. R. Lin, G. P. Guo, T. Tu, F. Y. Zhu, and G. C. Guo, Phys. 

Rev. Lett. 101, 230501 (2008); P. Pei, C. Li, J.S. Jin, and H.S. 

Song, e-print arXiv: quant-ph/101 1.2252vl (2010). 
[41] G. Burkard and A. Imamoglu, Phys. Rev. B 74, 041307 (2006). 
[42] J. M. Taylor and M. D. Lukin, e-print arXiv: |cond-mat/0605144| 

(2006). 

[43] J. M. Taylor, J. R. Petta, A. C. Johnson, A. Yacoby, C. M. Mar- 
cus, and M. D. Lukin, Phys. Rev. B 76, 035315 (2007). 

[44] O. N. Jouravlev and Y. V. Nazarov, Phys. Rev. Lett. 96, 176804 
(2006). 

[45] F. H. L. Koppens, C. Buizert, I. T. Vink, K. C. Nowack, T. Meu- 
nier, L. P. Kouwenhoven, L. M. K. Vandersypen, J. Appl. Phys. 
101,081706 (2007). 

[46] T. Yu and J. H. Eberly, Quant. Inf. Comp. 7, 459 (2007); M. 
Ali, G. Alber, and A. R. P. Rau, J. Phys. B: At. Mol. Opt. Phys. 
42, 025501 (2009); A.R.P. Rau, J. Phys. A: Math. Theor. 42, 
412002 (2009). 

[47] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998); S. Hill and 
W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 

[48] J. B. Yuan, L. M. Kuang, and J. Q. Liao, J. Phys. B: At. Mol. 
Opt. Phys. 43, 165503 (2010). 



